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1 - DEVELOPMENT O F MEASURING MEAN S 

1.1. - MEASUREMEN TS OF SURFACE FRICTION USING HOT WIRE 
GAUGES (Tests at F1/CFM) 


Six gauges were installed within the upper surface of an LC 
90 D airfoil, near the trailing edge. The flow is two -dim -■n- 
sional, velocity is constant and equal to to m/sec, the total 
pressure is equal to 1 bar and tne angle of attack vari^: 
between -5 and +13 degrees. Boundary layer profiles were 
measured for each configuration and at different ooints using 
the hot wire anemometer technique. The form parameters 
measured ' in that manner were used in the Ludwieg- Tillman 
formula in order to obtain the changes in the friction 
coefficient against which the values obtained from the gauges 
are compared. This comparison reveals good agreement between 
the Cf changes obtained by the two methods; the relative 
differences are in the order of 10‘/o, exceot in the 
neighborhood of the separation (indicated by an arrow) where 
the measurement technique is faulty. 




1.2. - MEASUREMENT OF THE SURFACE FRICTION, THREE-WIRE GAUGE 

To respond to the problem created by the measurement of the 

friction in unsteady boundary layers, various modifications 

were already added while building the wire gauges used in 

steady flow. These improvements make a large bandpass 

possible (much greater than 1000 Hz without an appreciable 

phase-shift) and also increase the sensitivity and the 

accuracy of the response of the Drobe in steady flow. The 

dispersion of the measurements perfo rmed does not exceed 3% 
*Numbers in margin indicate foreign pagination 
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of the friction. 


The determination of the sign of the friction can be obtained 
simply by adding two cold wires on each side of the sensitive 
element ; these wires make it possible to locate the warm 
wake . The measurement principle is the following: the rod 
. ovides tne friction coefficient ar.d the two cold wires give 
the temperature upstream and downstream from the rod. The 
unree signals being digitized, coT.rar.ison of the temperature 
data against a threshold provides the sign of the friction, 
ihe three- wire gauge • has been used to determine the 
distribution of t ^ over a cylinuer Dlaeed DerDendicula y 
to the wind. Flow fluctuations induced by the release of 
Karman vortices are clearly revealed by the surface friction 
measurements: Deriodic displacement of the stagnation ooint 

and a sign change of the friction in phase 0 = 180 degrees. 
On this last Doint, the instantaneous friction reaches rather 
large values, in the order of 3 to 5 Pa, whereas, for a reason 
of symmetry, the average value of must be equal to zero 
which is clearly shown by the signal processing that is 
performed . 

1.3. - MULTIPLE PROBES 

Hot-wire anemometer probes, rcultiwire probes (3 or 4 wire 
probes) have been used in order to define quickly and 
simultaneously the average and turbulent velocity fields of 
three-dimensional flows. 

Adapting these new measurement techniques to the industrial 
environment has been attempted in the FI wind tunnel of Le 
Fauga during tests to define vortex areas close to a sweptback 
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wing. Some difficulties inherent to the geometry of the probe 
itself and to the highly vortex-filled nature of the flow 
(rate of turbulence close to 50 %) have appeared and have led 
to the increase of the stiffness and strength of the sensor. 
On the otner ‘and , a new method for calibration and operation 
has been considered so as to oush back the limitations in 
tne use of the probe in turbulent flows tnat nave been highly 
deflected. 
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PREPARATION OF TESTS AT WIND TUNNEL F2 AT LE FAUGA 


This wind tunnel, used by GS£, is designed for the basic 
research work of the ONERA and DERAT Directorate of Aero- 
dynamics: complex three-dimensional configurations, the 

process of turbulence generation, the check of computation 
methods . 

i'he primary characteristics of the wind tunnel are: 

- atmospheric pressure, 

- velocity: 0 to 100 m/s 

- test section: 

w i d t h : 1 . 4 0 m 

height : 1.60 m 

length: 5 m 

- cross-section reduction ratio 12, low amount of 
tur oulence . 

The side walls of tne test section are easily removed. They 
are opaque or transparent depending on the need. The movable 
supDort of the laser velocity-measuring device that surrounds 
the test section provides a longitudinal displacement of 3.75 
meters, a vertical displacement of 1 meter and a transverse 
displacement of 0.6 meter from the aiming point, thus defining 
a large volume for test purposes. 

The first experiments will be performed by the DERAT and will 
have as parallel goals the refinement' of existing measuring 
systems (LDA in particular) and the definition of other means 
of measurement (hot wire exploration, for example). 

The assembly used will be a cylinder with a diameter of 12 
cm, motor-driven in rotation about its axis. It will make 
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it possible to conduct a fine and detailed study of the flow 
about the cylinder (average and fluctuating pressure, surface 
friction) with an accuracy of 0.1 degree by using a single 
generatrix of the cylinder. In particular, we will be 
interested in the area near the separation point and the 
stagnation point. Visualizations using a laser plane are 
also planned. 

A second thrust of the studies will be the cylinder wake. 
The size of the test section will indeed make it possible 
to explore the distant' wake (up to 40 diameters) and the 
characteristics of the wind tunnel will make it possible to 
conduct tests up to Renolds numbers of 700,000, that is to 
say beyond the critical Reynolds number. 

There again, visualizations using a laser plane are planned. 


3 - UNSTEADY BOUNDARY LAYERS 


COMPUTATION USING SMALL 


PERTURBATIONS 

For the computation of the incompressible two-dimensional 
unsteady flow boundary layers, a method using small 
oertur bat ions can prove very interesting by the rapidity witn 
which it can be apolied in all cases where the boundary 
conditions of the problem are periodic. Tne success of tnis 
type of method rests primarily on the fact that the non-linear 
effects are small, as has been shown exoerimentally or through 
integration of the complete boundary layer equations. 

The application of such a method in the more general case 
of a compressible three-dimensional boundary layer was 
necessary for certain practical applications of interest to 
industrialists such as the prediction of buffeting or flutter. 
Principle of the method 

To the three force, longitudinal momentum and transverse 
momentum integral equations, there are three primary unknowns 
H, 9 , 9^ which are the form parameters and the thicknesses 
of longitudinal and transverse momentum respectively. 


The system may be put in the form 


r. / EX • . 

I ( lj 


i!L + 


j = 1 


h^3x 


u e 3t 


) = E s 


SI 


i = 1, 2, 3 designates the equation considered and y. , j 

J 

= 1, 2, 3 represents the major unknowns, h is the length 
component in the x direction making an angle a with the 
outside velocity. The E coefficients are implicit functions 



Of U , a , 

e 


3ct 


hJx 


.9 , H, 0^ by means of a set of closing 
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relationships, and of the derivatives of ue in the z direction, 
making an angle x with respect to x. 

The small perturbation hypothesis consists of assuming that 
any f value can be written in the fcwm: 


f + f e 1 wt + f 


cc. 


-iiot with f the complex amplitude and 
c 


c 1 . 
<< 


f 

The system of equations then transforms itself into two 
systems of three equations, one for t: ; average values of 
the unknowns and not revealing fluctuating quantities, and 
the other for the complex amplitudes of the unknowns and a 
function of the average values. If the integration of the 
average system does not create a problem (it is identical 
to the system for the steady-flow problem), the integration 
of the unsteady system requires the computation of the complex 
portion of the E coefficients. The latter being a function 


of the independent variables, we have: E 


a E x. 


where 


- ? “ie ’ 

1 oi 

x . represents the complex amplitude of x. . The computation 
ie ■ 1 

of the partial derivatives of E relative to the x. 's can be 

l 

done algebraically in part or in a p.urely numerical manner 
by computing i values of E for increases of the xi's. The 
two procedures provide identical results but the second is 
faster in terms of machine time. Practically, the small- 
perturbation method requires only three times the amount of 
time for the steady-flow method. For one node in the mesh, 
0.053 second is required. 

In 2-D or 3-D incompressible flow the non-linearity effects 


remain very low and the method applies for A Ue/Ue = 0.20. 
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In incompressible flow (Me ^ 1.6), these effects become 

appreciable -as soon as AUe/Ue reaches several percentage 
points. 
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4 - BOUNDARY LAYERS OF BODIES OF REVOLUTION AND OF FUSELAGES 


4.1. - COMPUTATION OF OUTSIDE STREAMLINES 

The computation of the three-dimensional boundary layer 10 

requires that conditions be given at the limits, that is to 
say the magnitude and i. section of the velocity of the 

r i o n - v iscous f lo w . 

If the latter is determined from a perfect fluid computation, 
the velocity is simply broken down along the axes of the 
reference frame in which the boundary layer is comouted . 

On the otner hand, if the non-viscous flow is only known from 
the distribution of t l, e surface pressures, it is necessary 

to compute the direction of the velocity. The latter is 

determined in two steps. The s '•eamline passing through a 
given point is first computed in reverse: the direction of 
the velocity at that point is adjusted through an a firing 
method so that tne streamline ends at the stagnation point. 

This direction is then used to complete the solution 

downstream. 

This procedure has been applied to Airbus -type fuselages 
around 'which the pressure distribution has been obtained by 
means of a perfect fluid computation. The comparison of the 
directions of the velocity is entirely satisfactory. 

4.2. - BOUNDARY LAYER COMPUTATION 

The program for the three-dimensional boundary layer that 
was developed at the DERAT has been adapted to the case of 
fuselages (problems associated with the periodicity of the 
geometry and to the definition of transverse data at the ends 
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of the z grid). 

The program has been used for the computation of the boundary 
layer that develops over the nose of the Airbus-type 
aircraft. The separation at the bottom of the windshield 
is characterized by a large deflectuon of streamlines along 
the wall with respect to outside streamlines. 
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5 - PROBLEMS ASSOCIATED WITH WING-FUSELAGE INTERACTIONS 


5.1. - NUMERICAL STUDY 

The case of the Shabaka experiments was treated in detail 
in order to study the oerfortnances of the various turbulence 
models. A symmetrical wing without sweep or angle of attack 
is olaced on the floor of the wind tunnel. The computation 
begins 19 cm downstream from the leading edge of the wing, 
using the initial experimental data which indicate the 
presence of ' strong vortex. The equations are made Darabolic 
in the direction of the primary flow. The main problem 
consists of modelling the Reynolds stresses: 

Vortex viscosity model 


This is the simplest of all. The u'^u'^ 's are of 


the form: 


ti 2 . , 3Ui ^ 3Uj . 

u i“j= 3 S ij -v t ( — +—1-) 

3xj 3xi 
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k represents the turbulence kinetic energy and is 
the turbulent viscosity obtained, by example, through a 
formula for the length of the mixing. This model gives good 
results . 

" Com p 1 et e tr a ns po r t_ecj u a_t ion s_ f o r_ the _R e yn old_s_s_t£ e ss_e_s 
It is the most advanced model since it solves six 
equations for the six components of the Reynolds tensor : 


Du ' u * : 
— 1 j 


rij 


Dij 


ri.ll 

J 


3 Jijk 


Dt 3x k 

This model is defined as "complete" in the sense that all 
of the derivatives of the three components of the average 
speed are retained. The Launder-Reece-Rodi modelling is 
adopted for the redistribution term 0ij. The solution of 


24 



a transport equation for the dissipation makes it possible 
to close the system. The results obtained resemble those 
given by the vortex viscosity with, however, a certain 
displacement of the center of the vortex. 

System (6.6) 

In the preceding model, we neglect the convection and 

diffusion terms but we keep all of the aUi 's. At each 

axj 

computation node, it is necessary to solve i linear system 
of six equations for the six stresses. With respect to the 
results, there are few differences from the solution of the 


complete equations. 
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In the preceding system only the derivatives a ^ and a ^ are 
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retained. The Reynolds stresses are then expressed by means 
of simple analytical formulas, but the computation is in poor 
agreement with the experiments. 

Thus it appears that the derivatives of the secondary 
velocities V and W play important functions in the modelling 
of the u'i u'j 's. The vortex viscosity model, combining 
accuracy and a small number of computations, seems to be the 
most interesting but it cannot compute the development of 
the vortices created by turbulence. The model with six 
"complete" equations (or the system (6.6) in cases that are 
slightly convective), although more awkward to handle, offers, 
on the other hand, a greater range of applications. 
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5.2. 


STUDY OF THE MIXING OF A FUSELAGE BOUNDARY LAYER AND 


OF THE WAKE FROM A WING 

The experimental setup consists of an ONERA D wing airfoil 
with a chord length of 200 mm placed perpendicularly on the 
floor of the test section. Sweep and incidence angles are 
equal torero. 

The goal of the study is to study the three-dimensional 
mixture of the boundary layer on the floor which simulates 
a fuselage and of the wake generated by the wing. 

A first set of experiments had the ourpose of defining the 
conditions away i rom the floor, that is to say the development 
of the boundary layer and of airfoil wake. The velocity 
measurements have been accomplished over a rather tight grid 
with the help of a hot wire anemometer. The static oressure 
profiles were also measured. 

Then, measurements of the wake-boundary layer mixture were 
performed with the help of a probe with five pressure tubes. 
The mapping of the secondary velocities and the lines of equal 
velocity v/ere thus defined. They show in particular the rapid 
change in the nature of the flow near the trailing edge. 
In particular, the junction of the trailing edge and of the 
floor behaves as a hole for the secondary velocities. 

5.3. - EXPERIMENTAL STUDY OF THE FLOW IN THE NEIGHBORHOOD 
OF THE CONNECTION BETWEEN A FUSELAGE AND A SWEPTBACK 

WING WITH AN ANGLE OF ATTACK 

The experimental study consists of a wing airfoil (ONERA D 
airfoil with a chord length of 200 mm) placed with a large 
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sweep angle (60 deg.) and at an angle of attack (15 deg.) 

A measurement of the field of average velocities was performed 
with the help of a five-hole level indicator probe set within 
the angle defined by the wing and the floor of the test 
section . dine planes 'ere measured at X = -50 mm, 50 mm, 
100 mm, 130 mm , 200 mm, 250 mm, 300 mm, 350 mm, 332 mm (X 

is measured from the leading edge). 

.Measurements made with the help of a hot w're anemometer prooe 

(probe with one straight wire) were undertaken in order to 

validate the results obtained by the five-hole level indicator 

probe (with regard to the magnitude of only the velocity 

2 

vector) and also to define the u' component of the Reynolds 
tensor before undertaking a complete definition (6 components 
of the Reynolds tensor) of the turbulent field with the help 
of a four-hot wire anemometer probe. 
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6 - MODELLING . JRBULENCE 

6.1. - EFFECTS OF ROTATION ON TURBULENCE 

A general rotation of the flow, applied to an homogeneous 15 
and initially isotropic turbulent field, lessens the rate 
of decrease of the turbulent energy. However, this effect 
was not taken into consideration in the standard turbulence 
models . 

It had already been demonstrated that a sirr.Dle modification 
of the EDQNM two-point closing model makes it possible t j 
show the effect of rotation on spectral energy transfers, 
i h is & D Q h M model has been validated at very low Reynolds 

numbers by comparison with direct simulations of Navier-Stokes 
equations and at a moderate Reynolds number by comparison 
with simulations of large structures. In this latter case, 
the sub-node model describing the effect of small scales on 
large structures is derived from the EDQNM model, 
ihe next step consisted of performing a large number of 
numerical simulations of turbulence in rotation for a broad 
range of energy spectra, of Reynolds numbers, and of rotation 
rates with the EDQNM model. The tr.end analysis of the 
dissipation rate of the kinetic energy e as a function of 
time has provided a new model of the transport equation for 
this dissipation rate. Comparison with an experiment, in 
this case the experiment of Wigeland and Nagib, has made it 
possible to check the validity of the model. 

One significant point of the study of the rotation is that 
the shear encountered, inside the boundary layers for exam.ple, 
breaks down into a planar deformation and a rotation. Taking 

F. 

I - 
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the effects of rotation into consideration in the transport 


equations models of turbulent stresses has led to the 
improvement in predicting, both with cases of low shear and 
of high shear. 

6.2. - SPECTRAL INTEGRA! METHOD 

The spectral integral method is a new approach to formulating 
transport equations for turbulent quantities. It relies on 
a hypothesis defining the form of the sprectrum of the kinetic 
energy of turbr ence. For the time being, this model is only 
used to impro e the prediction of the dissipation rates of 
the kinetic energy e - 

In the case of a mono-scale model, the spectral integral model 
gives a new transport equation for e without an adjustable 
constant, contrary to the standard model. The trend of the 
dissipation rate is then tied to the production /d iss ipat ion 
balance and to the evolution of the very large structures 
that can be computed directly with a linear model. The 
prediction of the cases of homogeneous turbulence subjected 
to a flow is further improved. 

In the case of a multi-scale model, the energy spectrum is 
split into two areas: one describes the large structures 
carrying energy, the other describing the small structures 
of the inertial and dissipative areas. The spectral integral 
method makes it possible to replace the two transfer equations 
for the energy flow from the large structures toward the small 
structures and for the dissipation rate with one simple 
hypothesis in the form of the spectrum and one equation for 
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the energy ■ flux directly derived from the sub-node models 
developed at the DERAT. The multi-scale model is sensitive 
to the division of the production term between large and small 
structures. In the cases of pure rotation (no production), 
it gives very good results. Results close to and even better 
than those obtained v; ' th the mono-seals model are given with 
very intuitive divisions of production. 


TURBULENCE HOMOGENEOUS IN ROTATION 



Honeycomb Grid Test section in rotation 


Corrections of the equation for e 

* 0) l/2q 2 

Cl) - — - 1 — 

€ 

o computation without correction 
A proposed model 
B Bardina model 
HL Hanjalic and Launder model 



32 





spectral integral method 
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7 - DEVELOPMENT OF COMPUTATION METHODS 

7.1. - SOLUTION OF THE NAVIER-STOKES EQUATIONS PARABOLIZED 
IN THE TWO-DIMENSIONAL WITH SIMULTANEOUS COMPUTATION 
OF p , u , v . 

The Navier-Stokes equations are solved for two-dimensional 
laminar flow problems by means of a method called semi- 
elliptical. We take into account transverse variations of 

the longitudinal pressure gradient but the diffusion terms 
in the primary direction of the pressure gradient are 
neglected. In addition, . the longitudinal component of the 
pressure gradient is approached through a formula that is 
decentralized so as to recover the elliptical effects due 
to the pressure. 

The numerical method used is intended to solve the system, 
formed by the equations of motion and the continuity equation 
so that the velocity and pressure fields are comDUted 
simultaneously. The simultaneous handling of the equations 
indeed assumes a primordial significance in obtaining fast 
convergence and makes it possible to overcome difficulties 
such as separation at the trailing edge of an airfoil. Strong 
coupling between the viscous and non-viscous areas of the 
flow is automatically obtained. The system to resolve is 
written in the matrix form AX = B where X is a column vector 
having as components all of the variables of the system at 
the point where it is computed. Matrix A is a matrix that 
is tri-diagonal in blocks that we seek to inverse. The 
problem is solved step by step following the increasing values 
along the abscissa but the presence of the pressure terms 
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downstream from the computation abscissa implies the iteration 
of the process until convergence is obtained. This method 
requires only the storage of the pressure at all nodes in 
the mesh, the determination of U and V at the point where 

the computation is made implying the knowledge of these 
variables only at the previous point of measurement. This 

consideration together with . th e f acc them convergence takes 
place rapidly provides this technique with a defii.L-.e 
/ an tape in comparison with other methods that are 

classically used to solve this type of problem. 

The procedure described below has been checked for two 
different problems: 

The first of these studies refers to the flow around the 

trailing edge of a flat plate. The results of the semi- 
elliptical method are in excellent agreement with those 
obtained by classical coupling methods and with those obtained 
by the solution of the complete Navier-Stokes equations. 

The flow around a descending step is the object of the second 
application. Here, the longitudinal diffusion terms have 
been retained in the neighborhood of the step, downstream 
from it. The results are close to those obtained by more 
refined methods of the Navier-Stokes equations but the 
computation times are much less. 
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Channel with 



7.2. - BOUNDARY LAYER-PERFECT FLUID SIMULTANEOUS COMPUTATION 

WITH IN ’.R ACTION. APPLICATION TO SEPARATION BULBS. 

2(J 

For the computation of the flow about an airfoil, it is often 

advantageous to separate the viscous effects from the 

non-viscous effects. The area close to the airfoil will be 

governed by the boundary layer equations, the rest of the 

area being governed by the perfect fluid equations. Then, we 

will have to take into accoun , an interaction between these 

two areas. 

An attractive solution method consists of solving the two 
problems simultaneously ( Veldmann ) . Formally, the problem 
appears in the form: 

6 1 - CL (Ue) =0 and 6 1 - FP (Ue) = 0 
where CL and FP represent the boundary layer and perfect fluid 

operations linking the flow velocity Ue and the displacement 
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thickness 6^. Let us note that because of the elliptical 
nature of the non-viscous problem, an iterative process must 
be applied. 

This method has been applied to the computation of separation 
bulbs with transition for which a method of the inverse type 
has been developed . 

The boundary layer is computed with the help of an integral 
method for the la:.. _nar, turbulent or transitional regimes. 
The perfect fluid is considered as the result of the 
suoerposition of a known non-per turbed state and of a 
correction which is determined by computing a Cauchy 
integral. In that manner, the coupling problem involves only 
one space variable. 

The boundary and perfect fluid equations as well as the 
interaction rule connecting them are solved simultaneously 
by means of a relaxation method. 

The results have been compared against those obtained with 
the standard inverse method: they are very close to those 
results but provide a very considerable reduction in 
computation time. 


# 
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8 - TRANSITION OF THE BOUNDARY LAYER IN THREE-DIMENSIONAL 
FLOW • 

8.1. TRANSITION CRITERIA AND STABILITY 
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of low values of e (0 <e<10° ), the most unstable direction 

e^is responsible for the transversal instability, or rather 

quasi-transverse since it differs from the transverse 

direction by several degrees. The transition criterion is 

then obtained by an empirical correlation between, on the 

one hand, the Reynolds number formed with the displacement 

thickness of the velocity streamline projected in tne 

direction e . and, on the other hand, the form factor of 

m i n 


the longitudinal airfoil 


and the rate of outside 


turbulence: R6 (e=e . ) = g H (Tu). The application of thes< 


1 


rr.i.n 


1 


criteria in various series of experiments, including that 
of Boltz-Kenyon and Allen, has given satisfactory results. 

8.2. - COMPUTATION OF THE TRANSITION REGION 

The abscissa of the beginning of the transition being known, 
it is now necessary to compute the transition region itself. 
A standard method, known under the name of intermittence 
method, is used. A model relying on an equilibrium of the 
laminar and turbulent properties of the boundary layer has 
first been established by computing two-dimensional boundary 
layers through solution of local equations. It was then 
extended to the three-dimensional case. In addition, a first 
approach to a model for the solution of the integral equations 
of the boundary layer in the two and three-dimensional has 
been studied; it is already providing encouraging results. 


8.3. - PARAMETRIC STUDY OF THE SWEEP ANGLE EFFECT 
The OAP 01 airfoil (glider wing airfoil) is studied at three 
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On the upper surface, 


angles of attack: 0°, 1° and 2°. 

the distributions of the outside velocity are of the "laminar" 
type, with steep negative pressure gradients. There results 
from it a large transverse flow, and providing a sweep angle 
brings the transition closer to the leading edge. On the 
other hand, on the lower surface, the distributions are 
"flatter". The transverse flow, much weaker, has little 
chances to lead to a premature transition. Nevertheless, 
this transition exists but it is due to a longitudinal 
instability. The complete computation of the laminar, 
transitional and turbulent regions provides the variations 
of the thickness of the momentum at the trailing edge, 9 p p, 
which makes it possible to estimate the total drag coefficient. 

Parameter measuring is performed as a function of the 
angles of attack and of sweep. In that manner, we enow the 
effect of the sweep angle at a given angle of attack and a 
given Reynolds number. 
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(. i ) upper suriauc * '* (2) lower surface 


Distributions of C (0AP01 airfoil) 








8.4. EXPERIMENTAL STUDIES OF THE TRANSITION ON A SWEPTBACK 


WING 

Basic study at low speed 

A first phase consisted of pursuing the systematic study of 24 
l.oe effect of a sweep and of a Reynolds numoer. The sweep 
angle varied from 40 to 60 degrees and the normal angle of 
attack from -3 to -6 degrees, angles of attack for which there 
is no peak overspeed on the upper surface near the leading 
edge. For the -3° the pressure gradient is still negative 

whereas for the -6°angle of attack there is a velocity 
maximum at X/C = 0.56. 

The position of the transition is determined by means of two 
methods: visualization through sublimation and Jones criterion 
(measurement of the stagnation pressure in the boundary layer 
as a function of x at a constant altitude). These two methods 
show that for certain configurations the variation in 
magnitude of the transition abscissa can reach 10 to 15% of 
the chord length. These variations may be due to the fact 
that the infinite wing assumption is not checked (wing span 
too small in comparison with the width of the test section) 
or may be due to flaws in the model (geometry, roughness). 

This phenomenon is accentuated by the existence of steady 
flow instability waves approximately aligned with the outside 
flow and which appear on visualizations as "streaks". The 
detailed study of this phenomenon is under way. It has 
already shown the existence of a laminar boundary layer with 
characteristics that undulate greatly in magnitude in the 
area of the streaks. Good agreement has been obtained between 


the experimental determination of the wavelength and its 
prediction through a computation of laminar stability. 

The application of two transition criteria described in 
Paragraph 8.1 makes it possible to predict in a satisfactory 
manner the location of the transition caused by the transverse 
instability . 

Probings of the boundary layer made downstream of the 
transition area have shown that the downstream "undulation" 
of the boundary layer disappears when the latter becomes 
turbulent. The experimental values of the integral 
thicknesses of the boundary layer and of the friction 
coefficient are compared against a complete computation of 
the boundary layer, from the leading edge to the trailing 
edge. This computation makes it possible, in particular, 
to describe the transition area thanks to the intermittence 
method Dreviously described. 


Settling Chamber Test Room 
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Transition over a swept back wing in transonic regime 

A few preliminary tests aimed at obtaining data pertaining 26 
to the location of the transition over a sweptback wing in 
a moderately compressible flow have been performed at the 
T2 wind tunnel, operating at ambient temperature. 

They involve the’ "variable sweep" wing of airfoil LC 100 D 
with a chord length of 100mm; the results pertain primarily 
to the 45 deg. sweep angle. 

The transition was determined on the upper surface by means 
of visual visualizations using an oil film and oerformed at 
various negative angles of attack, with the upstream Mach 
number capable of reaching 0.3. 

We show the observed location of the transition as a function 
of the chord Mach number for various angles of attack, as 
well as the corres ponding results obtained by means of a 
computation of the transition using the transition criterion 
from transverse instability of incompressible flow. 

Despite the uncertainty occurring in the determination of 
the transition through visualization at the surface, the 
motion of the transition matches the one indicated by the 
criteria; we clearly note the forward motion of the 

transition when the angle of attack becomes more negative, 
a forward motion caused by the accentuation of the negative 
pressure gradients and the corresponding increase of the 
boundary layer transverse flow and its stabilizing effect. 
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CONDITIONS FOR THE APPEARANCE OF THE VORTEX REGIME 


OVER A SWEPTBACK WING 

Tests over the sweptback wing used in the basic study of the 
transition at a low speed have been devoted to both the study 
of the bursting of the leading edge transition bulb and the 
appearance of the vortex regime at high angles of attack. 

The passing from a flow with a leading edge bulb to a vortex 
regime has been revealed through the measurement of the 
pressures along the contour and through visualizations along 
the surface during a gr, ual reduction of the Reynolds number 
at a fixed angle of attack. 

The extent of the area affected by the bulb and the vortex 
flow area as a function of the Reynolds number could be 
determined . 

It is to be noted that artificially tripping the transition, 
performed during tests at a rather low Reynolds number, 
resulted in the suppression of the vortex regime. 
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PRESSURE COEFFICIENT DISTRIBUTIONS DOMAINS WHERE THE VORTEX 

AND BULB EXIST 
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9 - DETAILED STUDY OF THE CAST 7 AIRFOIL AT THE T2 WIND TUNNEL 


This experimental study dealt with the detailed analysis of 
the viscous flow regimes (boundary layer and near wake) near 
the trailing edge of the CAST 7 airfoil with a chord length 
of 200 mm, mounted between adjustable walls inside the test 
section of the T2 wind tunnel. 

All of the tests were performed with a transition tripped 
at 7 % of the chord, with a total pressure close to 1.7 bars, 
at ambient temperature. 

Two main configurations were considered for the tests: 

- a Mach number varying between 0.75 and 0.79 around the 
divergence Mach number at a zero incidence angle 

- an angle of attack varying between 1 and 3 degrees at 
Mq= 0.7, a configuration that causes the appearance of 

a local separation of the boundary layer at the bottom of 
a shock located between 20 and 40;4 of the chord. 

The measurements made inside the boundary layer were made 
primarily near the trailing edge, over the lower surface and 
particularly on the upper surface of the airfoil, a few 
explorations defining only the conditions upstream of the 
shock . 

The evolution of the displacement thickness and of the 
momentum thickness is given as an example of one case of the 
first configuration; we note a very quick thickening over 
the rear portion of the airfoil where the boundary layer 
develops into a velocity profile approaching a separation 
profile . 

A boundary layer computation using the 'experimental pressure 

t * 

r 
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distribution has been performed; there is good agreement with 
the experiment up to 75$ of the chord; a divergence then 
appears, particularly in the displacement thickness, and has 
led to a discussion on the possible effect of secondary 
factors (airfoil curvature, boundary layers on the side walls) 
to take into consideration in computations. However, we must 
emphasize the difficulty of accurately measuring a boundary 
layer close to separation near the trailing edge (vibration 
of the probe, inaccuracy in locating the surface). 

An example of stagnation pressure distributions measured at 
various Mach numbers in the near wake is also given; it 
clearly reveals the two areas of o^essurs drop originating, 
respectively, from the boundary layer and the shock wave, 
the second increasing clearly with the Mach number , that 
is to say in fact, with the strength and width of the shock 


wave . 



0 
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Displacement thickness and momentum of the boundary layer 
on the airfoil. 


55 




10 - SHOCK-BOUNDARY LAYER INTERACTION WITH TRANSITION AT 


THE BOTTOM OF THE SHOCK 

This experimental study dealt with the analysis of the 30 
behavior of the boundary layer and its . interaction with the 
recompression shock wave in the case where the boundary layer, 
initially laminar, becomes turbulent through transition inside 
the interaction region. 

The tests were performed at the T2 wind tunnel, operating 
at an ambient temperature, using the model of the CAST 7 
airfoil with a chord length of 200 mm. 

A first step consisted of defining the cases to be studied, 
searching for configurations providing a laminar boundary 
layer up to the shock; it involves pressure distributions 
that have, at their beginning, an area with a negative 
gradient stabilizing the boundary layer followed by an area 
with a zero or slightly negative gradient up to the shock. 

We chose to set the shock wave at about 50% of the chord; by 
changing the angle of attack and the Mach number, we obtain, 
ahead of the shock, a flow with a Mach number varying from 
1.18 to 1.40, thereby studying different strength values of 
the shock wave. 

The distributions -of the local Mach number show a two-step 
recompression with . a constant pressure plateau between the 
two steps which are characteristic of a lambda shock wave 
and of the existence of a boundary layer at the beginning 
of the interaction; the second recompression step, much 
clearer and much larger, is on the other hand characteristic 
of a turbulent boundary layer. 
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Surface visualizations using an oil film have in fact shown 
that the laminar boundary layer separates, transitions and 
returns to the surface in a turbulent flow near the inter- 
action . 

Measurements performed inside the boundary layer at the 
abscissa value X = 0 . 32 have confirmed that the boundary layer 
is laminar upstream of the shock; its characteristics , 
particularly its displacement thickness , depend 1 tie on 
the magnitude of the local Mach number. 

For each of the configurations, the boundary layer after the 
shock and its development as far as the trailing edge h a v ^ 
been described by means of measurements performed at abscissa 
va.lues of 0.65, 0.30 and 0.95. 

It is therefore clear that the boundary layer is turbulent 
after the shock and that the thickening that it nas undergone 
while crossing is greater with a higher upstream local Mach 
number. The displacement thickness then increases 
considerably down to the trailing ed~e, and even more so when 
the case considered involves a stronger shock. 

These experimental results provide the basis for checking 
a computation method that remains to ’be established and in 
'which the effects of a transition of the boundary layer on 
the boundary-shock wave interaction phenomenon are taken into 
account. The dashed line plotted on the figure only 
represents a very rough computed estimate, using the 
experimental distribution of the pressures with the hypothesis 
of a transition at the location of the shock; relatively close 
to the experiment with the weakest shock, the results of this 
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1 1 - ADJUSTING THE WALLS IN A TRANSONIC WIND TUNNEL 
11.1. - ADJUSTABLE WALLS FOR THE CRYOGENIC OPERATION OF T2 
In view of the cryogenic modification of T2 , we were led to 32 
equip the wind tunnel test section with new deformable walls 
operating at a low temperature. 

We chose to use thermally conducting metal sheets; thermal 
flow computations snow that it i. possible to have their 

temperature equal that of the flow over periods of time that 
are entirely compatible witn the duration of the blast, at 

least for a small thickness, and provided the effect of the 

straightened va.i_s is reduced to a pin - .at. 

a ! jus to ale alls of T2 are thus mace of Invar, a 
material with a cm-. aoa fflcieni. of ex ns ion , mechanically 
s v,u a -.0 i temperature; they have a thickness of 3 mm; 

tne straigntener supports are notched, with holes of a size 
computed to reduce their thermal inertia to the maximum 
extent. 

11.2. - OPTIMIZATION OF TECHNIQUES FOR ADJUSTING THE WALLS 
Four improvements were added to the original method: 

- Optimization and extension of the computation mesh to 

predict with greater accuracy the outside field using a 

formula by Green. 

- Determination of the velocity at infinity by cancelling 
the velocity correction, the weighing procedure giving 
greater weight near the model which, in fact, defines the 
Mach number of the test. 

Introduction of a projection making it possible to go 
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through the curved wall to make measurements along the 

straight reference line along which the outside field is 

* 

computed using Green's formula. 

- Extrapolation of the test zone up to x = + 40 m, by means 
of -a reference model (with distributed points) with a 

" connection *to the four outside experimental values. 

New numerical methods (finite elements, small transonic 
perturbations) have also been applied to compute the outside 
field in cases where the local Mach number at the wall is 
close to 1 . 

An analysis of the diagram of a velocity field through the 
use of a pannel method currently constitutes a first 
theoretical approach for the modification for three-dimen- 
sional flow. 

11.3. - EFFECT OF THE SIDE WALL BOUNDARY LAYERS 
A few tests of the CAST 7 airfoil between adjustable walls 
have revealed these effects in two ways: 

- a suction at a variable rate has made it possible to modify 
the thickness of the boundary layers on the upstream side 

* walls and their evolution inside the pressure gradients 
near the model; 

- the map of the displacement thickness around the airfoil, 
determined from . computations performed over different 
streamlines, has made it possible to build a detailed 
compensation over the cross-section, thanks to which the 
effects of the boundary layers seem to be greatly resorbed. 
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(3) Compensated side 
boundary layers 
Cn ) Suction ■ 

(5) Surface boundary layer 



Effect of the lateral boundary layers on the airfoil. 
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12 


DESCRIPTION OF THE CRYOGENIC FLOW INSIDE T2 


12.1. - THERMAL GRADIENTS INSIDE THE TEST SECTION 
A now analysis of the temperature distribution in space inside 
the test section was performed, a detailed measurement using 
a sh-ielded thermal probe examining the region close to the 
Wall in particular. 

The transverse distribution shows a uniform central area 
bordered by warmer flow areas that are clearly greater than 
the dynamic boundary layer; their size is practically 
independent of the test temperature and the difference between 
the wall temperature and . the stagnation temperature at the 
center does not exceed about 12 degrees Centigrade. 


12.2. - TEMPERATURE AND PRESSURE FLUCTUATIONS 


(Tests at M = 0.8, Pi = 2 bars, start of a throat) 

Pre_s£ure__ £l^ctuationsj_ They are measured up to a frequency 
in the order of 10 kHz by means of a Kulite sensor mounted 
on the wall of the test section (reduced frequency: n = 


fh 

U 


L 


n max 


F (n) dn = 

p (| u 2 ) 


r2 \ 2 ' 


At the ambient temperature, the spectrum is very close to 
the one that we had obtained before the circuit was modified. 
At a low temperature,, we note a slight rise of the fluctuation 
1 I which Lm certainly linked to the thermal turbulence. 


Temperature _ fluctuations They are measured inside the 
settling chamber, ud to a frequency in the order of 50 Hz, 
by means of a "cold wire". The average quadratic value of 
the fluctuation goes from 0.03K at ambient temperature to 
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O.lAK at Ti = 100K. Involving low frequencies, part of the 
fluctuation may come from the temperature control which is 
performed in steps; the homogeneity of the mixture of warm 
driving air and the driven cold gas may also be involved. 


T2 . 3 • - APPEARANCE OF PARTICLES - GAS CONDENSATION LIMIT 


Using the 

diffu 

sion 

of 

the light emitted by a 

laser beam, 

the thresholds 

when 

the 

particles appear as a 

function of 

temperature 

and 

the 

gas 

condensation limit were 

determined 


(tests at Pi - 2 bars). 

The recording of. the signal from the photomultiplier and the 
measurement of the pressure drop in the filter at the entrance 
to the settling chamber make it nossible to distinguish three 
phenomena : 

- appearance of ice particles at a temperature in the order 
of 250K then apnearance of carbonic snow particles at a 
temperature in the order of 1 3 5 K ; these particles "trapped" 
by the filter give rise to a greater pressure drop; 

- very large flow of particles at low temperature (Ti = 93K); 

the pressure and temperature conditions inside the test 

section make it possible to attribute it to the condensation 

of the air-nitrogen mixture; there appears that a total 

temperature of 100K can be easily obtained for transonic 

* / 

tests up to a total pressure reaching 3 bars. 
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Examnle of transverse temperature 


(distribution inside the test 


Stagnation temperature 
measurements near the 


section during a test at Ti=120K,; side wall. 


M=0.8 and Pi=2bars. 



Deser lotion of the pressure distri 
but ions at the test section wall . 


Evolution of the overall level 
of thermal fluctuation over 
the 0.5H 7 , to 50Hz frequency 
band as a function of the 
test temperature 
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Thu: figure shows the appearance of particles during an air 

blast at a low temperature Ti=95K. 
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13 - PREPARATION OF AIRFOIL CRYOGENIC TESTING AT T2 


-- 13.1. - SEQUENCE OF EVENTS FO R A CRYOGENIC BLAST 

The solution retained to reach a thermal equilibrium consists 
of precooling the airfoil before inserting it into the test 
secti-on. A cryogenic blast will thus include several 
'Sequences which will already have been perfected during tests 
conducted at T2 , without a model for the time being. 

- during the first steD, a flow at a low Mach number and at 

■ a low pressure but at a nominal temperature is initiated: 

the wind tunnel being started at the ambient temperature, 
a flow of liquid nitrogen at three times the equilibrium 
value is injected, causing the temperature to droD quickly; 
it is then gradually reduced, the computer providing control 
to reach the desired temperature level. 

- the second step consists of introducing the precooled 
airfoil, into the test section. 


- during the third step, the Mach number and the pressure 
are increased up to the values of the test itself through 
programmed action on the. flow of the air entering and on 

* the air exhaust rate. The temperature is maintained at 
a constant level through an increase of the nitrogen flow 
rate matching that of the flow of the air entering. 

- the fourth step, maintains the test parameters at nominal 
values thanks to a closed-loon pressure and temperature 
control; it is .during this step, where the conditions have 
been reached, that the measurements will be made on the 
airfoil. 
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The design ‘ and construction of a model of the CAST 7 airfoil 
for cryogenic tests at T 2 resulted from studies performed 
by the Direction des Grands Moyens d'Essais (Directorate of 
Major Testing Means) of ONERA, studies bearing on the 
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sy ste 

:m for me 

>ving 

the 

ore 

-cooling 

unit . 

A 

stub on 

the 

left 

sid 

:e makes 

it 

possible 

to lock 

the model 

when it 


is placed inside the test section. 

The chord length of the airfoil is 150 mm; it is equipped with 
103 pressure taps, the majority with a diameter of 0.3mm; in 
order to solve problems associated with the possibility of 
tripping the transition, 32 taps with a smaller diameter 
(0.8mm) are installed in the area of the leading edge where 
the boundary layer is thin. 

Temperature control is provided by 18 copper-constantan ther- 
mocouples placed at a depth of 1mm under the airfoil surface. 
The ten thermocouples are distributed along the center chord 
from the upper to the lower surfaces; the others are installed 
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along the wingspan. 


13.3. - PRECOOLING THE AIRFOIL * 

The preliminary phase of the test of a wing airfoil in a 38 
cryogenic flow therefore consists of cooling the model to 
bring it ta a temperature sufficiently close to that of the 
flow inside which it will be introduced. 

The solution retained is that of cooling by using a cold gas 
flow; a simple set-up having been bench-tested, it was used 
as a foundation for the Drecooling apparatus built for 
testing airfoils in a crvogenic regime in the T2 wind tunnel. 

This apparatus consists of a closed circuit operating at an 
atmospheric ambient ores sure and adjacent to one of the side 
windows in the test section; it includes two parts: 

- a cooling section hosting the airfoil and equiDped with 
the moving device which will make it possible to introduce 
it inside the test section of the wind tunnel, 

- the circuit itself, equipned w ith the nitrogen-injection 
system (two small-diameter, nozzles), the driving fan, the 

. exhaust, and the temperature control through a heater. 

The apparatus was first nerfected in the presence of a 
reference model made of an aluminum alloy; a satisfactory 

check of the time periods needed to reach different 

% 

tern, porn turn level r and cooling uniformity was made. 

It was then used to performed the precooling tests of the 
CAST 7airfoil Dlanned for the cryogenic tests at the T2 wind 
tunned . 

The examples of recordings given below show that the center 
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section of the model (region A) reaches the temperature of 
the flov; in about ten minutes; only the thermocouples located 
at 20 mm from the side st-ubs (region B) show differences 
caused by the effects of conduction, differences that should 
then disappear for the most nart during the cryogenic test 
itself. 


We were 

also able 

to check the 

good 

behavior of this model 

which i 

s of a new 

design , 

being 

made 

of a 

steel with a high 

n ickel 

content and 

cooled 

for 

the 

first 

time during these 


tests. No appreciable deformation via s recorded despite the 
application of temperatures dropping as low as 120K. Of 
courfie, proper behavior remains to be confirmed when the model 
is also subjected to aerodynamic stresses during testing. 


« 
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Legend : 

( 1 ) injection of 
airfoil 

(2) opening of 
digital valve 

; (%) 

(3) analog valve 
pre-set at 
38 % 



0 20 iOl 50 ao 




150 

100 

50 
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32 tans 0.1mm dia. 71 taps 0.3mm dia. 18 thermocouples 
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14 - PROBLEMS ASSOCIATED WITH THE USE OF CRYOGENIC WIND 
TUNNELS 

'14.1 - INSTRUMENTATION AND MEASUREMENT TECHNIQUES 

Pressure and temperature probes: In addition to a synthesis 40 
summarizing the results acquired with the various pressure 
sensors that’ we were led to test in a cryogenic environment, 
efforts were made to define and use miniaturized probes; a 
small-size probe (diameter of 4mm) and with a larger bandpass 
was built. 

A shielded stagnation temperature probe, with a recovery 
factor close to one, was also built; its utilization inside 
the wind tunnel T2 was very valuable when measuring the 
distribution of temperatures inside the wall boundary layer 
in the test section. 

Optical detection of oartiri.es: Preliminary tests inside the 
wind tunnel T'3 have shown the Dossibility of detecting the 
presence of particles inside the flow by using the light 
emitted by a low-power laser and diffused by the particles; 
the ontical viewing canability using two side windows that 
Veduce the formation of frost was necessary for "long" tests 
at T'3. The utilization of the system in the T2 wind tunnel 
has made it nossible to define the appearance of the different 
tynes of narticles % and to define the gas condensation limit 
in thn t, tunnel . 

Surface friction gauges: The DERAT hot wire and hot film 

gauges were designed and tested for mesuring the surface 
friction in low-temperature flows. A few tests conducted 
at T2 on a reference model (flat aluminum plate) have shown 
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the 

rod 

test 

the 

it 

and 


feasibility of 
gauge, capable 
ed inside the 
mounting onto 
to. define, in oar 
of the transition 


such surface measurements. The quartz 
of providing unsteady flow values was 
enclosure;* after a few modifications to 
the support, we expect to be able to use 
ticular, the nature of the boundary layer 
in the cryogenic mode. 


14.2. - T'3 CRYOGENIC WIND TUNNEL WITH FAN 

The improvement of the operation of the small T’3 cryogenic 
facility was pursued with the installation of a new exhaust 
casing now located right after the test section outlet; the 
suction of the boundary layer at the entrance to the diffuser 
contributes in a major fashion to a reduction in the pressure 
drop of the circuit and to the expansion of the operating 
domain of the wind tunnel. 

The preparation of a two-dimensional set-up particularly 
designed for the study of the detection of the transition 
in a cryogenic environment was performed; this set-up 
includes a sim.Dlified model composed of a flat plate with 
.an unsymmetr ical leading edge and an adjustable flap, and 
computed to minimize oversneed at the seam. 


The 

test section is equipped 

with simplified 

deformable walls 

that 

make it possible 

to 

avoid blockage 

and to adjust the 

P f’ 

ure distribution 

a 1 on p; 

the flat pla 

te. A system for 

t 

longitudinal and 

vertical exploration 

of the flow is 

planned, particularly 

for 

flat plate 

boundary layer 


measu remen ts . 
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Detailed view of side windows 
Legend : 
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(4) Compressed air 

(5) Side walls 

( 6 ) 

(7) Lenses 

(8) Photomultiplier 
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( 1 0 ) F i 1 1 e r 

! (ll)Output (voltmeter) 
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j High voltage 


Terr.nerature probe. 
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Optical detection of particles 
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Legend : 

- ( 1 ) Low velocity exhaust 

. (2) Hin;h velocity exhaust 

( 3 ) M. , 

test .section 

(4) Test section 

(5) Throat 

(6) Liquid nitrogen injector 



T'3 wind tunnel 
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